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ABSTRACT 

cn ■ 

o : 

^ ' We study the cluster of H I and O VI absorption systems and the claimed 

,_H ■ detection of O VIII absorption from the intergalactic medium at z ~ 0.0567, 

1^ ■ associated with a group of galaxies toward the BL Lac object PKS 2155-304. 

' As measured by spectrographs on the Hubble Space Telescope, Far Ultraviolet 

^ ■ Spectroscopic Explorer, and Chandra, this system appears to contain gas at a 

CN ■ variety of temperatures. We analyze this multi-phase gas in a clumpy-infall 

' model. From the absence of C IV and Si III absorption in the Lya clouds, 

O ■ we infer metallicities less than 2.5-10% of solar values. The only metals are 
O ■ 

■ detected in two O VI absorption components, offset by ±400 km s from the 
^ ' group barycenter [cz ~ 16,600 km s~^). The O VI components may signify 
O ■ "nearside" and "backside" infall into the group potential well, which coincides 

• with the claimed O VIII absorption. If the claimed O VIII detection is real, 
our analysis suggests that clusters of strong Lya and O VI absorbers, associated 

2 ■ with groups of galaxies, may be the "signposts" of shock-heated, metal-enriched 

^ ■ baryons. Through combined UV and X-ray spectra of H I and O VI, VII, and 

■ VIII, one may be able to clarify the heating mechanism of this multiphase gas. 

■ Subject headings: intergalactic medium — quasars: absorption lines — ultravio- 
■ — ' let: general 
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1. INTRODUCTION 

With the development of high-throughput uhraviolet and X-ray spectrographs, we can 
begin to account for a significant fraction of the low-redshift "missing baryons" (Fukugita, 
Hogan, & Peebles 1998; Shull 2003) in the intergalactic medium (IGM). Recent surveys by 
the Hubble Space Telescope (HST) and the Far Ultraviolet Spectroscopic Explorer (FUSE) 
show that 32 ± 6% of the baryons reside in the warm (10^ K) Lyct forest (Penton, Stocke, & 
Shull 2003), while 5 - 10% he in the shock-heated (10^-^ K) IGM traced by O VI absorption 
(Tripp, Savage, & Jenkins 2000). A predicted 30-40% of the IGM remains to be found in 
even hotter gas (10^~^ K) through X-ray absorption lines of O VII, O VIII, and perhaps 
Ne IX (Fang et al. 2002; Nicastro et al. 2002). 

The total baryon density is known to 10% accuracy from measurements of hght-element 
nucleosynthesis (D/H) and acoustic oscillations in the cosmic microwave background (CMB). 
Recent estimates from these techniques yield consistent values for Qh, the fractional contri- 
bution of baryons to the closure density: Qb = (0.020 ± 0.002)/i-2 (D/H - Buries & Tytler 
1998) and Qb = (0.0224 ± 0.0009)/i-2 (CMB - Spcrgcl et al. 2003; Ncttcrficld ct al. 2002) for 
a Hubble constant Hq = lOO/i km s^^ Mpc^^. Even Hq has become a well-measured param- 
eter, with values oi h — 0.71 ± 0.08 (Freedman et al. 2001) from the Cepheid Key Project 
with HST. Using the CMB values, the co-moving baryon density is pb = (4.2 ± 0.2) x 10~^^ 
g cm~^ and the hydrogen density is tih — (1.90 ± 0.2) x 10~^ cm~^ for primordial hehum 
abundance Yp — 0.244 by mass. 

Numerical simulations of the low- 2; IGM (Cen & Ostriker 1999; Dave et al. 2001) predict 
that the gas is distributed nearly equally (30 ± 10% each) in three "phases": (1) warm 
photoionized gas (10^ K); (2) warm/hot shocked gas (10^~^ K); and (3) collapsed halos, 
galaxies, and clusters with much hotter gas (T > 10^ K). The warm phase is observable 
through UV absorption lines of H I, C HI, C IV, Si HI, and Si IV, while the warm/hot 
phase is detectable in O VI (UV) and in X-ray lines of O VII, O VIII, and Nc IX. The 
warm-hot intergalactic medium (WHIM) is thought to be produced by shock-heating during 
gravitational infall into dark-matter filaments co-evolving with the rise in cosmic metallicity. 
These heavy elements are theorized to be expelled from galaxies by tidal stripping (Gnedin 
1998) and "starburst winds" as observed by Martin et al. (2002). The metals are then 
assimilated into infalling clumps of IGM. Details of this model remain unclear. Are the 
heavy elements expelled from galaxies in hot or warm gas? How are the metals mixed with 
the IGM? What is the extent of their transport? Addressing these issues requires detecting 
and correlating many warm and hot IGM absorbers. 

In this paper, we provide clear evidence for multi-phase IGM toward the BL Lac object 
PKS 2155-304. This gas is associated with a small group of galaxies at z ^ 0.0567 and 
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a cluster of 7 strong Lya absorbers (ShuU et al. 1998) at velocities cz — 16, 283 — 17, 570 
km s~^. This environment exhibits the diffuse warm phase (Lya, Ly/3, Ly7), the shocked 
(WHIM) phase (O VI), and unconfirmed evidence (O VIII) for much hotter, virialized gas 
residing at the bottom of the group's gravitational potential well. 

In § 2, we describe our FUSE and HST observations of these absorbers in H I and 
O VI (and limits on C IV, Si III) and compare them to the claimed O VIII absorption. In 
§ 3 we examine the absorber kinematics for the possibility of nearside/backside infall and 
ionization-state (O VI/0 VIII) consistency. We conclude in § 4 with a discussion of an infall 
model that may be consistent with large-scale structure in this group. In this picture, the 
H I and O VI absorbers arise from clumps of gas falling at 200-400 km s^^ relative to a much 
hotter (lO^-^-'^' o K) substrate that could be detectable in O VII or O VIII absorption. The 
metaUicities may range from 0.02-0.1 solar. All the phases predicted by simulations may be 
present in this group, suggesting that clusters of strong Lya and O VI absorbers associated 
with groups of galaxies may be the "signposts" of shock-heated, metal-enriched baryons. 



2. OBSERVATIONS 

In HST Cycle 9, we obtained new data from the STIS Echelle (E140M) with 10 orbits 
(PID 8125, 28.5 ksec) and compared the absorber velocities to those in a small group of 
galaxies containing four large H I galaxies and one dwarf seen with the VLA (see Fig. 3 in 
Shull ct al. 1998). The PKS 2155-304 sightline passes near the center of this group, none of 
these galaxies lies closer than 400 kpc to the sightline, and no other dwarf galaxies were seen 
down to iriB ~ 19. The galaxy ID velocity dispersion is (Xgai = 325 km s^^, and the Lya 
velocity dispersion is CLya = 740 km s~^. We estimate the galaxy overdensity as 5 ^ 100, 
based on the group of five galaxies. 

We also observed PKS 2155-304 with FUSE (Moos et al. 2000) in three visits on 23-24 
October 1999 (P1080701, tcxp = 19.2 ksec, and P1080705, tcxp = 38.6 ksec) and 18 June 
2001 (P1080703, tcxp = 65.4 ksec). All data were calibrated with CALFUSE 2.2.1. We base 
most of our analysis on the second visit, which obtained data for all of the FUSE detectors, 
whereas the earlier visits occurred before the SiC and LiF channels (Sahnow et al. 2000) were 
aligned. For the LiFla segment (1000-1084 A), we roughly tripled the amount of available 
data by co-adding night-only reductions of the three observations. This method recovered the 
redshifted Ly7 absorption (1027.6 and 1028.0 A) fr om below O I airglow emission (1025.76 
and 1026.47 A). We applied small wavelength offsets (< 0.05 A) to the FUSE data to align 
narrow interstellar lines and place the FUSE and STIS data on a common velocity scale set 
by Galactic 21- cm emission. 
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The claimed X-ray detection of O VIII absorption (Fang et al. 2002) with 4.5(7 signifi- 
cance was made with the Chandra Low-Energy Transmission Grating Spectrometer. How- 
ever, the redshifted O VIII absorber is not seen near 20 A in the spectrum taken by the 
X-Ray Multi-Mirror Observatory (XMM) (see Fig. 1 and Table 1 in Rasmussen, Kahn, & 
Paerels 2003). Although the XMM velocity resolution is poorer than that of Chandra /YiETG, 
the XMM data have high S/N and detect oxygen absorption at 2; = with equivalent widths 
16.3 ± 3.3 mA (O VIII, 5.7 a significance) and 9.0 ± 2.7 mA (O VII, 8 a significance). The 
non-confirmation of the z — 0.0567 O VIII absorption remains a puzzle. 

Figure 1 shows a velocity overlay of the Ly«, Ly7, O VI, and O Vlll absorbers between 
cz = 15,000 — 20,000 km s^^. The detected lines and their properties are summarized in 
Table 1; the three strong Lya absorbers are labeled A, B, and C for clarity. We detect Ly/3 
at 16,243 km s'^ (A), but Ly/3 hues at 16,973 km s"^ (B) and 17,110 km s"^ (C) are severely 
blended with the N II A1084 interstellar hne. We detect Ly7 in components B and C, but 
we derive a large error on their equivalent width due to continuum placement uncertainty. 
We also detect O VI in components A and C, at 16,243 and 17,144 km s~^. 



3. INTERPRETATION OF THE SPECTRA 

With our new HST/STIS E140M data, we can derive more accurate column densities 
and metaUicities for the strong Lya absorbers. In our previous study (ShuU et al. 1998) at 
20 km s"^ (HST/GHRS) resolution, we estimated Nri = (3-10) x 10^^ cm'^ for the blended 
B and C components at 2; ~ 0.0567. There, we set a metallicity limit Z < 0.003^© from the 
absence of Si 111 and C IV absorption. At the higher S/N and better spectral resolution of the 
ST1S/E140M, together with FUSE measurements of Ly7, the derived H I columns decreased 
considerably. The integrated Lya optical depths give N(H 1) = 1.0, 2.2, and 0.94 x lO^'' 
cm~^, respectively, for absorbers at 16,185, 16,974, and 17,147 km s~^, consistent with the 
observed Ly7 strengths. 

We still did not detect any Si 111 A1206 or C IV A1548 absorption at the Lya velocities 
(Table 1). The 6-12 mA equivalent width limits correspond to N(C IV) < 2.9 x 10^^ cm^^ and 
N(Si III) < 2.7x 10^^ cm~^. As described in ShuU et al. (1998), we assume that the absorption 
arises in a uniform slab of gas photoionized by the metagalactic radiation background, with 
spectral shape for a background dominated by AGN with intrinsic spectral index of a = 1.8 
modified by IGM absorption (Shull et al. 1999). The ionization correction (conversion from 
C IV/H 1 to [C/H]) depends only on the ionization parameter, logU oc Iq/hh, where Jq is the 
specific intensity at 1 ryd and Uh is the volume density of the gas. For —2.0 < logf/ < —0.5, 
the metallicity upper limit from C IV/H I ranges between 0.025-0.10 of the solar abundance 
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in the CLOUDY photoionization code, (C/H)0 = 3.55 x 10"^ (Grevesse & Noels 1993). A full 
discussion of the parameter dependence of the metallicity estimate appears in a subsequent 
paper. 

The two detected O VI systems have column densities of N(0 VI) = (2.8 ± 0.5) x 10^^ 
cm~^ (A) and (3.3 ± 0.5) x 10^^ cm~^ (C) assuming a linear curve of growth. At the 
temperature (log Tqyi = 5.45) of maximum ionization abundance (&ovi = 17 km s~^), the 
stronger O VI A 1032 lines in both components would still be unsaturated. We identify these 
two O VI absorbers with Lya systems at 16,185 (A) and 17,147 km s^^ (C). No O VI is 
seen in the 16,974 km s~^ (B) absorber. The O VI lines are the only metals yet detected 
in these absorbers, but their metallicity is uncertain, owing to an indeterminate ionization 
state and thermal phase (see Savage et al. 2002). If the O VI arises in photoionized gas, its 
column density and hmits on C IV imply log U > —0.8, for C/0 in a solar abundance ratio 
of 0.50 ± 0.07 (AUende Prieto, Lambert, & Asplund 2002). 

If the ionizing background has Iq ^ 10~^^ ergs cm~^ s~^ Hz~^ sr~^ (ShuU et al. 1999), the 
density constraint is uh < 10'^-^ cm~^, and the corresponding hne-of-sight extent of the gas is 
D > 1.1 Mpc, comparable to the extent of the galaxy group. We beheve that this low-density, 
large-absorber model is highly unhkely. Multiple absorbers of this extent could not fit within 
the 1 Mpc volume of the group, without undergoing collisions and shredding. However, if 
the O VI arises in collisionally ionized gas, the requirement that the gas be so rarefied no 
longer holds. In addition, the fraction of O VI in collisional ionization equilibrium peaks at 
T ~ 10^'^ K, and it is tempting to associate the O VI with the WHIM. The tentative detection 
of O VIII absorption makes collisionally ionized O VI still more compelling. The claimed 
strength of the X-ray feature (Fang et al. 2002) requires log N(0 VIII) 16.0 ± 0.2, which 
could represent hot gas at the group barycenter. However, redshifted O VIII absorption was 
not confirmed by XMM, nor has O VII absorption at 2; = 0.0567 been detected (Nicastro et 
al. 2002). 



4. DISCUSSION 

In order to understand the cluster of PKS 2155-304 absorbers ai z ^ 0.0567, one must 
resort to a multiphase model. The H I, O VI, and O VIII systems are unhkely to co-exist, 
owing to kinematic offsets and physical considerations. For example, the claimed O VIII 

absorption centroid (16, 624 ±237 km s~^) appears to differ from the observed H I and O VI 
absorbers (A, B, C components). Also, the observed Ly7 line widths (FWHM < 30 km s~^) 
rule out H 1 temperatures greater than 20,000 K. The measured O VI line widths (40 ± 10 
km s~^ FWHM) make it unlikely, although marginally possible within errors, for O VI and 
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O VIII to exist at the same temperature, logT 6.25 ±0.1, where (O VI/0 VIII) ^6x 10'^ 
in colUsional ionization equihbrium (Sutherland & Dopita 1993). However, we discount this 
scenario, because of the observed velocity offsets (~ 400 km s^^). For similar reasons, the 
strong H I absorbers cannot co-exist with the hot O Vlll. At the high temperatures needed 
to ionize O VIII coUisionally, the H I ionization fraction would be too small to be detected. 

Therefore, we favor a more complex scenario in which the H I and O VI absorbers arise 
in clumps of gas falling into a small-group potential. Any hot gas, visible in O Vll or O VIII, 
would then exist at the barycenter of the group (16,600 km s^^) which would be consistent 
with the kinematic offsets. If this model is correct, the two O VI absorbers (A and C) 
would arise from clumps undergoing "backside" and "nearside" infall, respectively. Because 
production of coUisionally ionized O VI requires shock velocities Vs > 130 km s"^ (ShuU 
& McKee 1979; Dopita & Sutherland 1996), the infall model has consequences for which 
Lya absorbers contain associated O VI. In particular, the absence of detectable O VI in the 
strongest Lya system (B) may be a result of a low relative velocity between the infalling H I 
clump and the substrate. From Figure 1, one sees that O VI is present in the H I absorber 
at 17,147 km s~^, located approximately 170 km s^^ redward. Thus, one might speculate 
that these two blended Lya absorbers, one with and one without O VI, represent shocked 
and unshocked gas. 

In the clump-infall model, the relative velocity between the two O VI components and 
the group barycenter is ~ 400 km s^^, which would prodiicc a post-shock temperature of 
~ 2 X 10^ K, similar to that (2.5 x 10^ K) at which the observed O VI and claimed O VIII 
would co-exist in coUisional ionization equilibrium. The cooling time of such gas is 



where n_4 = n///(10~^ cm~^) and Tg.s = (T/10^'^ K). Here, A_23 is the radiative cooling rate 
coefficient in units of 10'^^ erg cm^ s^^, typical of gas at 10^'^ K with 0.1 solar metallicity 
(Sutherland & Dopita 1993). In the range 5.0 < logT < 7.0, A(T) cx T^^, so tcooi is longer 
than the age of the universe, unless the gas temperature can be decreased (Tg.s ^ 1) or 
the density increased (n_4 ^ 1). Through shock-compression and cooling, the cooling will 
accelerate as A(T) moves toward its 10^-^ K peak. 

We can relate the hydrogen density to the cosmological overdensity 5 by: 



(20 Gyr)r6.5 

n_4A-23 



(1) 



UHiz) = (1.90 X 10-^ 



0.046 



{1 + 6). 



(2) 



At the redshift {za — 0.0567) of the strong absorbers, we have n_4 (5/500). If the infalling 
gas is shocked to T, = (1.24 x 10^ K){Vs/300 km s"^)^ to form the WHIM phase, it must 
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acquire a fairly high density {S > 500 at z ^ 0.1) in order to cool, recombine, and form the 
observed O VI. Adiabatic shocks can provide a compressive trigger, but further radiative 
cooling is required, since adiabatic compression lengthens the cooling time (tcooi n^^^ for 
T oc n^/^ and A(T) oc T~^). Lower- velocity shocks, Vg ~ 200 km s~^, would cool more 
rapidly, with temperatures near the peak of the cooling curve. 

Two additional issues concern the ratios of C IV, O VI, and H I. In the two detected 
O VI absorbers, (C IV/0 VI) < 10~^'^, whereas models of radiatively cooling gas (Dopita 
& Sutherland 1996; Indebetouw & Shull 2003) find (C IV/0 VI) ^ 0.1 - 0.2. Where is 
the C IV? Is the cooled, rccombined layer truncated in the hot-gas environment? Similarly, 
we do not understand the wide variation in [N(H I)/N(0 VI)] ~ 0.1-10 among other IGM 
absorbers (Shull 2003). The production of shocked O VI requires relative velocities above 
~ 150 km s~^ and 0/H abundances above a few percent solar. However, the production 
mechanism of O VI remains unclear (Sembach et al. 2003). While shock-heating is the 
preferred mechanism, other processes have been proposed involving interfaces between hot 
substrates and moving clouds - thermal conduction, shear instabilities, and turbulent mixing 
(see Indebetouw & Shull 2003 for a review). These processes all have problems explaining 
the observed components, N(0 VI) ^ 3 x 10^^ cm~^. Each conductive interface typically 
produces N(0 VI) ~ 10^^ cm~^, while shock heating requires supersonic flow and depends on 
the Mach number cubed. If Thot ~ 10^'^ K (adiabatic sound speed Cs ~ 240 km/s), infalling 
clumps at 300-400 km s"^ have Mach numbers of just 1.25 - 1.67. 

Resolving the O VI heating paradox may require understanding the cloud-substrate 
interfaces at a deeper level. Alternatively, if the O VIII absorption turns out to be non- 
existent, a virialized hot substrate may not have formed yet. Inter group gas at T < 10^ K 
would produce sufficient shock-heating, and the observed O VI absorption may arise from 
the interactions of infalling clumps with metals in a intergroup medium injected from the 
galaxies. With combined UV and X-ray spectra, one may be able to distinguish among these 
processes through the amount of shock heating and the relative abundances of O VI, O VII, 
and O VIII. 

This work is based on data obtained for the Guaranteed Time Team by the NASA- 
CNES-CSA FUSE mission operated by the Johns Hopkins University. Financial support to 
U.S. participants has been provided by NASA contract NAS5-32985. Wc were also supported 
by grants GO-08571.01-A from the Space Telescope Science Institute and NAG5-7262 from 
NASA/LTSA. We thank Taotao Fang for providing the Chandra spectrum that appears in 
Figure 1 and John Stocke and Bill Blair for helpful discussions. 
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Table 1. PKS 2155-304 IGM Absorption Lines 



Line 


veiocixy 




L\ V 




( Kill S j 


(lUAj 


( Kill S J 


H I Lya (A) 


16185 d= 5 


331 ± 12 


75 ±8 


H I Lya (B) 


16974 ± 5 


495 ± 10 


76 ±5 


H I Lya (C) 


17147 ±5 


353 ±8 


110 ± 10 


H I LyP (A) 


16243 ± 10 


85 ± 10 


92 ±17 


H I Ly7 (B) 


16936 ± 10 


16 ± 4 


< 30 


H I Ly7 (C) 


17109 ±10 


13 ± 3 


< 30 


H I Ly6^ 




<12 


< 30 


VI 1032 (A) 


16243 ± 10 


37 ± 7 


35± 10 


VI 1038 (A) 


16252 ± 10 


25 ± 5 


30± 10 


VI 1032 (C) 


17144 ± 10 


44 ± 6 


40± 10 


VI 1038 (C) 


17116 ±10 


Fe II blend 




VIII 18.969 


16624 ± 237 


u.otli 


< 1380 


C III 977^ 




< 10 




C IV 1548"^ 




< 12 




Si III 1206'^ 




< 6 





'^Full width at half maximum (FWHM) from 
STIS/FUSE profile fitting, or by converting Chandra line 
width of a{0 VIII) < 0.039 A (Fang et al. 2002). The 
doppler parameter, 6dopp = FWHM/ [2 (In 2)^/2]. 

^Upper limit at 4cr significance. 
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Fig. 1.— Absorption lines of H I (Lya from HST/STIS, Ly7 from FUSE), O VI 1032/1038 
(FUSE), and claimed detection of O VIII with Chandra (Fang et al. 2002) toward PKS 2155- 
304. The three strongest Lya absorbers are labeled (A,B,C - see Table 1). We detect O VI 
in the 16,243 km s~^ (A) and 17,144 km s~^ (C) absorbers, and Ly7 in B and C absorbers; 
claimed O VIII absorption is offset by ~ 400 km s~^ from Lya and O VI (vertical dotted 
lines). Interstellar lines (Ly/3, C II, O VI) are labeled in panel 3. Lower curve in panel 
4 shows O VI 1038 offset from O VI 1032. The Fe II line shows an asymmetric profile 
that suggests presence of O VI 1038 at 17,140 km s~^, but blending precludes an accurate 
measurement. 



